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Abstract. We present new high angular resolution observations at near-IR wavelengths of the core of the Luminous 
Blue Variable n Carinae, using NAOS-CONICA at the VLT and VINCI at the VLT Interferometer (VLTI). The 
latter observations provide spatial information on a scale of 5 milli-arcsec or ~11 AU at the distance of rj Carinae. 
The present-day stellar wind of r\ Carinae is resolved on a scale of several stellar radii. Assuming spherical 
symmetry, we find a mass loss rate of 1.6xl0 -3 M0/yr and a wind clumping factor of 0.26. The VLTI data taken 
at a baseline of 24 meter show that the object is elongated with a de-projected axis ratio of approximately 1.5; the 
major axis is aligned with that of the large bi-polar nebula that was ejected in the 19th century. The most likely 
explanation for this observation is a counter-intuitive model in which stellar rotation near the critical velocity 
causes enhanced mass loss along the rotation axis. This results from the large temperature difference between 
pole and equator in rapidly rotating stars. r\ Carinae must rotate in excess of 90 per cent of its critical velocity 
to account for the observed shape. The large outburst may have been shaped in a similar way. Our observations 
provide strong support for the existence of a theoretically predicted rotational instability, known as the Q limit. 



Key words. Stars: circumstellar matter: stellar winds: 
mass loss, stars: individual: 77 Carinae 

1. Introduction 

The Luminous Blue Variable ij Carinae is the most lumi- 
nous star known in the galaxy l|Davidson fc Humrihrevsl 
119971) . Its extreme properties make it an interesting lab- 
oratory to study the physics of the most massive stars 
in galaxies. Not much is known about the life of such 
stars, including their birth and post-main-sequence evo- 
lution, n Carinae has already left the main sequence and 
is now in an unstable phase; it experienced a large out- 
Send offprint requests to: R. van Boekel: 
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* Based on observations obtained at the European Southern 
Observatory 



burst in the 19th century which created the beautiful bi- 
polar nebula seen in many images (referred to as the ho- 
munculus). As m uch as 10 may have been ejected 
during that event IjSmith et al.l2003a) . A smaller outburst 
seems to have occurred around 1890, creating a smaller, 
but similarly shaped nebula hidden inside the larger one 
l|lshibashi et al.ll2003|) . 

The origin of the highly bi-polar shape of the homuncu- 
lus is a matter of debate. Model calculations show that a 
spherical explosion into a mas sive equatorial toru s can ex- 
plain the observed geometry l)Frank et al.l H995L Indeed, 
evidence for the pr esence of a 15 t orus was found from 
ISO spectroscopy l|Morris et al.ll999|) . Note however, that 
ISmith et al.l ( 2003b) argue that most of this mass is ac- 
tually located in the lobes. Other models reproduce the 
shape of the nebula by assuming a non-spherical outburst 
that runs into a spherical envelope. Recently, several au- 



2 



Van Boekel et al.: Direct measurement of size and shape of r\ Carinae wind 



N 




— 0.8 L 

0.6 0.4 0.2 0.0 -0.2-0.4-0.6 
A RA (arcseconds) 



Fig. 1. Diffraction- limited 2.29 fim image of the core of the 
r\ Carinae nebula, taken with the adaptive optics camera 
NAOS-CONICA at the VLT. The position of the Weigelt 
Blobs is indicated. 

thors proposed that luminous stars rotating close to their 
critical speed have stellar winds with a higher wind d ensity 
and expansion velocit y at t he poles (Owocki et al . Il996l 
Maeder & Desiacaues l200lL Dwarkadas & Owocki 12002ft . 
Therefore, the shape of the homunculus may be a natural 
consequence of the combined effects of rapid rotation and 
the high luminosity of 77 Carinae. Note that the extreme lu- 
minosity of 77 Carinae implies that even a modest amount 
of rotation cau ses the star to be close to its critical velocity 
l|LanEferlll99"7|) . Evidence in support of a polar enhanced 
wind was recently inferred from Hubble Space Telescope 
(HST) spectroscopy of starlight scatte red off dust grains 
in the homunculus (Smith et al. l2003ah . 

In this Letter we present the first results of an extensive 
high angular resolution near-IR study of the core of the 
homunculus, revealing for the first time the shape of the 
present-day stellar wind on a scale of 5 milli-arcsec. 

2. Observations 

We have observed 77 Carinae with the adap tive optics 
camera NAOS-CONICA (Lenzen et al. Il998l Rousset et 
al. Eb00) attached to Yepun, one of the 8.2 meter Unit 
Telescopes of the Very Large Telescope (VLT) of the 
European Southern Observatory (ESO), located at Cerro 
Paranal, Chile. NAOS was in the visual wavefront sen- 
sor configuration with 14 x 14 subapertures used for 
wavefront sensing. CONICA was used with the S13 cam- 
era (13.25 milli-arcsec/pixel) and the NB_239 intermedi- 
ate band filter with a central wavelength of 2.39 /im and a 
width of 60 nm. We employed a neutral density filter with 
an attenuation factor of «70 at the observing wavelength 
in order to avoid saturating the central peak of the point 



spread function. Individual exposures were 0.16 s long, and 
20 exposures were co-added, resulting in a total exposure 
time of 3.6 s per frame. Four such frames were obtained. 
The reduced image, diffraction limited at 70 milli-arcsec, 
is shown in figure Q 

The central 1.5 arcsecond region of the nebula is dom- 
inated by a point source and shows a complex mor- 
phology of blobs in the immediate vicinity. The bright- 
est blobs were previously discovered u sing s peckle imag- 
ing tec hnique s (Weigelt & Ebersberger 1986, Hofmann & 
Weigelt 1 1988|). There is a clear asymmetry in the emission, 
which is brightest in the north-west. Approximately 57 per 
cent of the flux seen inside a region with a diameter of 1.4 
arcsec, centered on the core of the nebula and correspond- 
ing to the Airy disk of the siderostats, is concentrated in 
the unresolved point s ource. We use recent literature data 
iSmith fc Oehrd EoOO) to flux calibrate our observations 
and find that the unresolved central source has a flux of 
about 200 Jy. 

We used the two 35 cm test telesc opes of the VLT 
Interferometer (G lindemann et al. 120031) and the test in- 
strument VINCI IjKervella et alJl2000l) to obtain intcrfcr- 
omctric measurements at baselines ranging from 8 to 62 
meters in length. The observations were carried out in the 
first half of 2002 in four different nights, and again in early 
2003. For a full description of the observations, and anal- 
ysis of the data we refer to Scholler et al. (submitted to 
A&A). The baselines used, have a ground length of 8, 16, 
24, and 66 m respectively. In particular observations with 
the 24 m baseline cover a wide range of projected baseline 
orientations. 

In figure we show the VINCI observations, given as 
squared visibility points. The visibility curves were cor- 
rected for the emission in the siderostat beams on scales 
larger than 70 milli-arcsec, by adding the resolved compo- 
nent seen in our NAOS-CONICA images. This results in a 
reduction of the visibilities by a factor (l-f res ), where f res is 
the fraction of the total flux in the siderostat beams that 
does not come from the central source (0.43). For each 
baseline, the top panel of figure shows the average visi- 
bilities at that baseline, while the bottom panel shows the 
variations of visibility with projected baseline orientation 
for the 24m baseline. 



3. Analysis 

3.1. Nature of the emission 

Figurc|21shows that the visibility decreases with increasing 
baseline, indicating that the interferometer is resolving the 
central source. The near-IR emission seen in the central 
region of 77 Carinae must come from an object that emits 
half of its flux in an area of 5 milli-arcsec or 11 AU diame- 
ter, assuming c ircular symmetry and adopting a distance 
of 2300 parsec IjDavidson fc Humphrevslll997P. Using th e 
measured size and flux of 200 ,Iv llSmith fc Gehrd fcOQOL 
a lower limit to the temperature of the object can be de- 
rived, assuming that at a wavelength of 2.2 /jm there is 
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no extinction. We find a minimum temperature of 2300 K. 
Allowin g for a foreground extinction of one (two) mag- 
nitudes 15mie7eiaL||200l|) , this temperature goes up to 
3200 (5000) K. The minimum temperature is too high to 
be caused by thermal emission from dust. We conclude 
that we have spatially resolved the ionized stellar wind. 

3.2. Mass loss rate 

We constructed a simple physical model for the observed 
size and flux of r\ Carinae, assuming that it is due to a 
spherical star with a dense, ionized and isothermal stellar 
wind that reaches a terminal flow velocity of 500 km/s 
(Hilli er et alJl200l|) . The wind is further assumed to be 
clumpy, with the clumping factor / defined as ~p = / p 
and where p is the unclumped wind density. For a more 
detailed description of the model assumptions, we refer to 
de Koter et al. (in preparation). 

A range of mass loss rates and clumping factors 
are consistent with the data. When combined with op- 
tical spectra of th e central star taken with the HST 
l|Hillier et a,lJl200lh . a unique determination of the mass 
loss rate and the wind clumping factor is possible. The 
HST spectrum can be fitted using a non-local thermo- 
dyna mic equilibrium (non-LTE) model for the star and 
wind (Hillier c tai]l200l|) which is degenerate in the ratio 
M/yff, where / is the wind clumping factor as defined 
above. A good fit to the HST observations is obtained for 
a ratio M /y/J = 0.00316 M /yr. We finally arrive at a 
mass loss rate of 1.6±0.3x 10~ 3 M Q /yr and a wind clump- 
ing factor / = 0.26 (figures [5] and |3J). The spatial light dis- 
tribution of the non-LTE model yields a somewhat better 
fit to the data than the LTE model. 



3.3. Asymmetry of the emission 

The 24 m baseline (figure |2J) data show a variation of visi- 
bility with projected baseline orientation, which is incon- 
sistent with a spherical distribution of light at 2.2 /zm. 
Deviations from spherical symmetry will introduce addi- 
tional uncertainty in the derivation of the mass loss rate 
given above. To determine the axis ratio of the emission, 
a gaussian shape was fit to the observations at each pro- 
jected baseline orientation. Note that the actual light dis- 
tribution is not gaussian. Therefore, the typical size of 
7 milli-arcsec derived from the gaussian fits differs from 
the 5 milli-arcsec derived from our (spherical) physical 
model. We find a major to minor axis ratio of 1.25 ± 
0.05, and a position angle of 134 ± 7 degrees. This po- 
sition angle is within errors equal to the position angle 
of 132 degrees l)Davidson et ahlboOlf) found for the ho- 
munculus. We therefore have conclusive evidence that the 
density contours in 77 Carinae's wind are elongated along 
the major axis of the Homunculus. The chance that the 
two prolate structures - i.e. the core and the Homunculus 
- have different orientation in 3 dimensions but still show 
the observed alignment when projected on the plane of 
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Fig. 2. upper panel: Observed visibility as a function of 
projected baseline. The full line is the best fit LTE wind 
model to the VLTI observations. The dashed line is the 
best fit non-LTE radiative transfer model to the observed 
visibilities. Lower panel: Variation of FWHM fitted to the 
visibility of rj Carinae measured with VINCI at the VLTI 
as a function of projected orientation of the 24 m base- 
line. The solid line gives the best fit to the measurements, 
assuming a 2D gaussian shape of the source at each pro- 
jected baseline orientation. The amplitude of the size vari- 
ations gives a ratio of major to minor axis of 1.25±0.05. 
The major axis has a position angle of 134±7 degrees east 
of north. 

the sky is about 10 percent, given the uncertainties in 
our measurements. The core on a scale of 10 1 AU, and 
the Homunculus reaching out several times 10 4 AU, must 
be truly aligned in space. Adopting an angle between the 
major a xis of the Homunculus and the line of sight of 41 
degrees l|Davidson et al.ll200l|) . the observed aspect ratio 
of major over minor axis of 1.25 implies a de-projected 
aspect ratio of roughly 1.5 (in detail this de-projection 
depends on the actual three dimensional structure of the 
object). 

4. Discussion 

It has been t hought (e.g. Lamers & Pauldrach ll99li Poe & 
Friend 1986) that stellar rotation enhances mass loss in the 
equatorial regions, resulting in disk-like winds. However, 
this would imply that r\ Carinae's rotation axis is per- 
pendicular to the major axis of the bi-polar homunculus, 
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Fig. 3. Allowed ranges in mass loss rate and wind clump- 
ing factor for the LTE wind model (contours) and the full 
non-LTE radiative trans fer model applied to the Hubble 
Space Telescope spectra (Hil lier et al.ll200l|) (dashed line). 
Both data sets agree for a mass loss rate of 1.6xl0 -3 
]VI©/y r and a wind clumping factor of 0.26. 



which is unlikely. A recen t model (Owocki et al. fl996l 
Dwarkadas & Owocki 2002) for line-driven winds from lu- 
minous hot stars rotating near their critical speed predicts 
a higher wind speed and density along the poles than in 
the equator. This counter- int uitive effect is cau sed by the 
increased polar temperature l|von Zeipe]lll924h and asso- 
ciated radiation pressure. Our data clearly favour the po- 
lar wind model. The VLT data do not provide informa- 
tion abo ut the velocity field . However, recent HST spec- 
troscopy l|Smith et all2003a|) of starlight reflected by dust 
in the Homunculus indicates a latitude-dependent wind 
outflow velocity, with the highest velocities near the pole; 
this is expected for a wind which is stronger at the poles. 
These data also suggest a polar enhanced win d den sity. 
Applying the model of D warkad as & Owocki l|2002L see 
also Maeder & Desiacaues l200l|) . we find that rj Carinae 
must rotate at about 90 per cent of its critical velocity to 
account for the observed shape. 

The question arises whether the model assump- 
tions (line-driven wind, radi ative envelope) made by e.g. 
iDwarkadas fc Owockil l|2002l) are applicable to rj Carinae. 
We note that y C arinae is almost certainly fully convec- 
tive l|Langerlll997|) . due to its near-Eddington luminosity 
(defined as the luminosity where surface gravity is com- 
pensated by radiation pressure). Therefore the difference 
between polar and e quatorial te mperatures in rj Carinae 
will likely be smaller l|Lucvll967|) than in polar wind mod- 
els, that adopt radiative envelopes. For such a convective 
envelope to produce a substantial temperature contrast, 



r\ Carinae must rotate in excess of 0.9 of the critical speed. 
The observed mass loss rate of 1.6xl0 -3 M /yr can be 
explained in terms of radiation driven wind theory (C. 
Acrts, private communication). 

The alignment of the homunculus and the present-day 
wind suggests a common physical cause. Rotation may 
then also be responsible for the shape of the homunculus 
(|Dwarkadas fc Owockll2002(L An outburst in 1890 prob- 
ably produ ced a bipolar nebula with a present-day size 
of 2 arcsec l|lshibashi et al.11200^ . which is aligned with, 
and inside the larger homunculus. There is thus strong 
evidence that the wind geometry is similar over a wide 
range of mass loss rates. It is not likely however that 
line-driven wind models are applicable to the outbursts. 
Nevertheless, our data underpin the importance of rota- 
tion for the post-main-sequence evolution of very massive 
stars such as r\ Carinae; it seems inevitable that as the star 
evolves, it will run into a rotationa l instability, referred to 
as the fi limit l Langer e t al. 1999). 



Acknowledgements. We would like to thank C. Aerts, P. Morris 
and N. Langer for useful discussions. DJH acknowledges partial 
support from NASA grant NAGW 3838. 

References 

Davidson, K. & Humphreys, R. M. 1997, ARA&A, 35, 1 
Davidson, K., Smith, N., Gull, T. R., Ishibashi, K., & Hillier, 

D. J. 2001, AJ, 121, 1569 
Dwarkadas, V. V. & Owocki, S. P. 2002, ApJ, 581, 1337 
Frank, A., Balick, B., & Davidson, K. 1995, ApJ, 441, L77 
Glindemann, A., Algomedo, J., Amestica, R., et al. 2003, in 

Interferometry for Optical Astronomy II. Edited by Wesley 

A. Traub. Proceedings of the SPIE, Volume 4838, pp. 89-100 

(2003)., 89-100 

Hillier, D. J., Davidson, K., Ishibashi, K., & Gull, T. 2001, 

ApJ, 553, 837 
Hofmann, K.-H. & Weigelt, G. 1988, A&A, 203, L21+ 
Ishibashi, K., Gull, T. R., Davidson, K., et al. 2003, AJ, 125, 

3222 

Kervella, P., Coude du Foresto, V., Glindemann, A., & 
Hofmann, R. 2000, in Proc. SPIE Vol. 4006, p. 31- 
42, Interferometry in Optical Astronomy, Pierre J. Lena; 
Andreas Quirrenbach; Eds., 31-42 

Lamers, H. J. G. & Pauldrach, A. W. A. 1991, A&A, 244, L5 

Langer, N. 1997, in ASP Conf. Ser. 120: Luminous Blue 
Variables: Massive Stars in Transition, 83 — h 

Langer, N., Garcfa-Segura, G., & Mac Low, M. 1999, ApJ, 520, 
L49 

Lenzen, R., Hofmann, R., Bizenberger, P., & Tusche, A. 1998, 
in Proc. SPIE Vol. 3354, p. 606-614, Infrared Astronomical 
Instrumentation, Albert M. Fowler; Ed., 606-614 
Lucy, L. B. 1967, Zeitschrift fur Astrophysics, 65, 89 
Maeder, A. & Desjacques, V. 2001, A&A, 372, L9 
Morris, P. W., Waters, L. B. F. M., Barlow, M. J., et al. 1999, 

Nature, 402, 502 
Owocki, S. P., Cranmer, S. R., & Gayley, K. G. 1996, ApJ, 
472, L115+ 

Poe, C. H. & Friend, D. B. 1986, ApJ, 311, 317 

Rousset, G., Lacombe, F., Puget, P., et al. 2000, in Proc. SPIE 

Vol. 4007, p. 72-81, Adaptive Optical Systems Technology, 

Peter L. Wizinowich; Ed., 72-81 



Van Boekel et al.: Direct measurement of size and shape of r\ Carinae wind 



5 



Smith, N., Davidson, K., Gull, T. R., Ishibashi, K., & Hillier, 

D. J. 2003a, ApJ, 586, 432 
Smith, N. & Gehrz, R. D. 2000, ApJ, 529, L99 
Smith, N., Gehrz, R. D., Hinz, P. M., et al. 2003b, AJ, 125, 

1458 

von Zeipel, H. 1924, MNRAS, 84, 665 

Weigelt, G. & Ebersberger, J. 1986, A&A, 163, L5 



